Anomalous Roughness in Dimer Type Surface Growth 
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We point out how geometric features affect the scaling properties of non-equilibrium dynamic 
processes, by a model for surface growth where particles can deposit and evaporate only in dimer 
form, but dissociate on the surface. Pinning valleys (hill tops) develop spontaneously and the surface 
facets for all growth (evaporation) biases. More intriguingly, the scaling properties of the rough one 
dimensional equilibrium surface are anomalous. Its width, W ~ L a , diverges with system size L, 

. This originates from a topological 



as a = i instead of the conventional universal value a = 
non-local evenness constraint on the surface configurations. 
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The theory of non-equilibrium dynamic statistical pro- 
cesses has developed rapidly in recent years. Driven sys- 
tems display intriguing scaling properties and can un- 
dergo various types of dynamic phase transitions |Q. 
Kardar-Parisi-Zhang type surface growth is an exam- 
ple There, the properties of the depositing (evap- 
orating) particles are not specified, but are implicitly 
presumed to be geometric featureless monomers. In sur- 
face catalysis type processes the geometric shape of the 
molecules matters. The onset of the catalytic process 
is associated with a so-called absorbing state dynamic 
phase transition ^-^|. Monomers give rise to directed 
percolation and dimers to directed Ising type transitions. 
Subtle geometric features are known to be important in 
equilibrium crystal surface phase transitions as well. The 
competition between surface roughening and surface re- 
construction depends on topological details of the crys- 
tal symmetry. Those determine whether a reconstructed 
rough phase can exist or not Geometric features are 
also important in diffusing particle systems. The shape of 
diffusing particles introduces conservations and leads to 
anomalous decay of particle density autocorrelations . 
Therefore, the natural question arises, whether and how 
the shape of the deposited particles influences the growth 
and equilibrium properties of crystal surfaces. 

Consider a crystal built from atoms of type X. As- 
sume that deposition always takes place in dimer form, 
X2, aligned with the surface. The dimer attaches to 2 
horizontal nearest neighbour surface sites. It looses its 
dimer character after becoming part of the crystal. Evap- 
oration can take place only in Xi molecular form, but a 
different partner is allowed. In this letter we study the 
one dimensional (ID) version of this process. This can 
apply to step shapes during step-flow type growth on 
vicinal surfaces. The adsorbed particles do not diffuse in 
this version of our model. However, topological features 
that drive our results are preserved, even when monomer 
diffusion is allowed but limited to terraces. Jumps across 



steps are unlikely due to Schwoebel barriers ||. So our 
main results do not alter in systems with diffusion. 

We describe the ID surface configurations in terms of 
integer height variables h{ — 0,±1,±2, They are 
subject to the so-called restricted solid-on-solid (RSOS) 
constraint, hi~hi + \ = 0, ±1, and periodic boundary con- 
ditions, /iL+i = hi. The dynamic rule is as follows. First, 
select at random a bond + I). If the two sites are 
not at the same height, no evaporation nor deposition 
takes place. If the two sites are at the same height, de- 
position of a dimer covering both sites is attempted with 
probability p, or evaporation of a dimer with probability 
q = 1 — p. Processes are rejected if they would result in 
a violation of the RSOS constraint. 

The first surprise is that the surface always facets 
during growth and evaporation, although the surface is 
rough in equilibrium. The second surprise is that the 
equilibrium surface width W ~ L a scales with an anoma- 
lous small exponent a = 0.29 ± 0.04. The data could be 
consistent with an even smaller value due to the strong 
finite size scaling corrections in Fig. 0(a). 

ID surfaces, irrespective of being in equilibrium or in 
a stationary growing (evaporating) state, display, with 
only a few very specific exceptions, the universal rough- 
ness exponent a = h. The up-down aspect of the steps 
become uncorrelated beyond a definite correlation length, 
and therefore the surface roughness obeys random walk 
statistics at large length scales, which implies a = \. 
Our dissociating dimer deposition process circumvents 
this universal argument by means of a novel type of non- 
local topological constraint. The dimer aspect implies 
that all surface height levels must be occupied by an even 
number of particles. However, due to the dissociative 
nature of the dimers that information is not preserved 
locally. The "evenness" constraint is non-local. At local 
length scales the surface looks the same as in monomer 
deposition processes, but the global surface is much less 
rough. We checked this numerically. Define a window of 
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length b. The surface roughness scales as W ~ b a , with 
a = i for b <C L, but crosses-over to the global finite size 
scaling exponent a ~ 0.29 ± 0.04 for b — ► L. 

We performed a detailed numerical study of the prop- 
erties of even- visiting random walks || that are globally 
restricted to visit each site an even number of times. The 
results, together with analytical scaling arguments, yield 
the value a = This value lies within the numerical 
error bars in Fig. 0(a) for the dimer deposition model. 

The details of our random walk study are rather tech- 
nical and will be presented elsewhere , but the essence 
can be captured by the following intuitive scaling argu- 
ment. Consider the even-visiting random walks for time 
interval < t < T(= L). We assign a defect variable to 
each site, to mark that it has been visited by the random 
walker an odd/even number of times up to time t. The 
even-visiting constraint is satisfied when all defects dis- 
appear at time T. Initially, the random walker does not 
feel the constraint and diffuses freely for t < Tf ree <C T. 
The defects are uniformly spread over a region of size 

1/2 

£ ~ r /ree ■ Then it stops spreading and starts to heal the 
defects. The healing time for a single defect in the region 
of size £ is order of £ 2 . By assuming that the healing 
process for each defect is independent, we estimate the 
total healing time T/ leQ ; ~ £ d + 2 with spatial dimensional- 
ity d. As T hea i > T/ ree , we conclude that £ ~ T l ^ 2+d \ 
i.e., a — 1/3 for d = 1. Existence of a time scale Tf ree 
explains the crossover behavior of the surface roughness 
in the window length b. Similarly, we can argue that the 
surface width scales with a = ^ in generalized n > 2 
dissociating n — mer type deposition processes. 

The dynamic critical exponent z at the equilibrium 
point follows from how the surface width diverges as 
function of time, W ~ tr . We find numerically that 
(3 = 0.111 ± 0.002, see Fig. 0(b). This suggests the value 
z = 3, since z — a/ (3. 

The equilibrium surface roughness is unstable with re- 
spect to growth and evaporation. It facets immediately. 
This phase transition is second order. The correlation 
lengths that characterize the faceted structure diverge. 
Before addressing this issue we need to describe and ex- 
plain the faceted phase. The valleys in the growing sur- 
face are sharp and the hill tops rounded, see Fig. |[ This 
shape is inverted for p < q. The faceting is caused by the 
spontaneous formation of pinning valleys during growth 
(pinning hill tops during erosion, for p < q). Consider 
for example, dimers on a flat surface for p > q. Odd seg- 
ments between them act as the nuclei of pinning valleys. 
Such valleys can not be filled by direct deposition. The 
only way to get rid of them is by lateral movement of the 
sub- hills. 

In finite systems the surface grows in shocks. An initial 
rough or flat configuration grows fast at first, but pinning 
valleys appear randomly at all surface heights. The inter- 
face develops into a rough faceted structure with many 



sub-hills and growth almost stops. From here on the sur- 
face advances only when sub-hills anneal out by the lat- 
eral movement of the pinning valleys. The annealing time 
of a sub-hill scales exponentially with its size. This expo- 
nentially slowing down healing process leads ultimately 
to a faceted W shape with only two remaining pinning 
valleys. Their lifetime diverges exponentially with the 
lattice size. After their demise the surface experiences a 
growth spurt, and the cycle restarts all over. 

The mechanism for lateral movement of pinning val- 
leys is exchange of active bonds between ramps. Active 
bonds are locations along the ramp where a dimer can de- 
posit or evaporate. Most steps on the ramps are only one 
or two atomic units wide and therefore dynamically inac- 
tive, see Fig. ^. Active bonds move up or down the ramps 
by deposition or evaporation of dimers. The growth bias 
p > q gives them an upward drift. This must lead to an 
exponential distribution. 

Fig. ||(a) shows the logarithm of the active bond dis- 
tribution, p{x), versus the horizontal distance x from the 
center of a hill top for various values of p > q. The 
straight lines for large x confirm the exponential distribu- 
tion of active bonds along the ramps, p(x) ~ exp[— 
We determined this from odd lattice sizes, in particular 
L = 257, where the surface contains an odd number of 
pinning valleys and therefore reaches a V shaped station- 
ary state in which it remains pinned at all times. 

Every now and then an active bond moves in the op- 
posite direction, against the flow, and reaches the valley 
bottom. That pinning valley moves by two lateral lat- 
tice constants when the active bond jumps across onto 
the other ramp. The probability for this is very small, 
and scales exponentially with the ramp length, but it is 
larger from the lower ramp, and therefore the valley bot- 
tom moves in the direction of the lower hill, and actually 
accelerates, because that lower hill keeps shrinking. 

Near the rounded hill tops, the surface remains highly 
active and initially the active bond density does not de- 
crease significantly with x. This defines a second charac- 
teristic length scale, £oi representing the flatness of the 
rounded hill tops, see Fig. 0. Surprisingly, both lengths, 
£/ and £o, diverge at p = q. Fig. ||(b) shows that the 
curves in Fig. ||(a) collapse according to a single length 
scaling form 

p{e,x) = e^T{e v x) (1) 

with v = 1.0(1) and p = 0.0(1). This means that on ap- 
proach of the p = q critical point the hills maintain their 
shape in the sense that £o and £f diverge simultaneously 
and with the same exponent £o ~ £/ ~ (p — q) • 

It is surprising that both length scales of the faceted 
phase diverge at the equilibrium point. The structure of 
the rough phase would be much more complex when one 
of them remained finite. This actually happens in the 
following generalization of the model. 
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Recently, Alon et. al. |T^] added to the conventional 
monomer type RSOS model growth the constraint that 
evaporation from flat segments is forbidden. It remains 
unclear how this can be experimentally implemented, but 
the interesting aspect of their model is the presence of 
a roughening transition, belonging to the directed per- 
colation (DP) universality class S, and unconventional 
roughness properties at this transition. Our model be- 
comes a directed Ising (DI) type generalization of 
this when we disable digging on flat surface segments. 
Modify the evaporation probability q to rq when both 
neighbours are at the same height as the update pair 
+ 1), i.e., hi-\ = hi = h i+ i = h l+2 . At r = 0, the 
no-digging limit, it becomes impossible to dig into the 
crystal layers beneath the current lowest exposed level. 
That level itself becomes frozen as well when it fills-up 
completely. 

Figure |J shows the phase diagram. The rough equi- 
librium point broadens into a rough phase (the shaded 
area). Along the DI ~ E phase boundary, the surface 
growth is zero. Inside the rough phase the surface grows, 
but slowly. Its scaling properties are complex and ob- 
scured numerically by huge corrections to scaling. The 
surface width W seems to grow logarithmically in time 
for L > 2 10 , and maybe the stationary state width scales 
logarithmically as well, but extracting the true scaling 
properties is a rather hopeless endeavor. 

The origin of this complexity is easy to pin point. The 
rough surface grows slowly, although the bare coupling 
constants p < q favours evaporation. It performs a del- 
icate balancing act. The surface erosion picture from 
r = 1 still holds along faceted ramp segments. There the 
surface evaporates due to the downward drift velocity of 
active bonds along the slope, but this is frustrated by 
the emergence of pinning hill tops. The surface grows 
at flat surface segments due to an upward pressure cre- 
ated by the reduced digging probability factor r, but the 
formation of pinning valleys limits this. Moreover, the 
non-local evenness constraint is at work as well. Growth 
and evaporation are dynamically balanced only along the 
faceting transition line DI — E (Fig. |J). Everywhere else 
the rough surface grows slowly. 

The properties of the two faceted phases confirm the 
above intuitive picture. The erosion faceted phase sig- 
nals the local stability of eroding ramps inside the rough 
phase, while the growth faceted phase indicates that flat 
segments persist. We have numerical evidence showing 
that the active bond characteristic length of the ero- 
sion faceted phase does not diverge along the roughening 
transition line DI — E for r < 1. This confirms that 
eroding ramps remain locally stable. On the other side 
of the phase diagram, the flatness length scale £o of the 
growth faceted phase does not diverge along the p = q 
roughening line. This confirms the persistence of locally 
stable flat segments inside the rough phase. 

To illustrate this, we present some of the details of 



the latter. Recall that along r = 1, the active bond 
distribution p(x) for different p collapses onto a single 
curve (Fig. ||(&)). For r < 1 this fails. At the transition 
point p = q, p(x) scales algebraically, as p(x) ~ 
see Fig. ||(a), but only beyond the central flat part of 
the hills. The flatness length scale, £o remains finite. Its 
value varies as £o ~ |1 — r \~ v ^ w i tn v — 1-0(1) along the 
line p = q. In Fig. ||(a), £o can be as large as £o — 40. 
This explains the poor finite size convergence of the sur- 
face roughness inside the rough phase. 

The scaling properties of the p = q faceting transition 
follow from the behaviour of the active bond density. Nu- 
merically, the surface width scales at p = q as W ~ L, 
like in the faceted phase. From the perspective of the 
rough phase the faceting transition takes place when the 
total active bond density, p — L p(x)dx, vanishes, be- 
cause in the faceted phase the £o segments of the rounded 
hill tops are of measure zero compared to the ramp seg- 
ments. We find p ~ (q — p)Pp with f3 p very close to 1. 
At p — q itself, the powerlaw p{x) ~ a; -1 predicts that 
p scales with system size as p ~ In L/L. The numeri- 
cal data in Fig. ||(6) confirm this. Finally, p decays at 
p = q algebraically in time with exponent 0.32(1). This 
suggests a dynamic exponent z = 3.1(2). 

Erosion below the currently lowest exposed level be- 
comes strictly forbidden at r = 0. The evaporation 
faceted phase becomes flat, and a directed Ising (DI) type 
roughening transition takes place at p = pui = 0.317(1). 
Hinrichsen and Odor already documented this. They 
independently introduced the r = limit of our model. 
They also report that the surface width scales at poi as 
y/\og(t), and inside the rough phase as log(f). 

In summary, the presence of a non-local topological 
constraint on equilibrium surface configurations in dis- 
sociating dimer type surface growth, leads to anomalous 
reduced surface roughness, with exponent a = 4 instead 
of the conventional value a = \. Moreover, the grow- 
ing (evaporating) surface is always faceted, due to the 
spontaneous creation of pinning valleys (hill tops). Un- 
der other circumstances, in particular when the digging 
probability on flat surface segments is being suppressed, 
an intermediate slowly growing rough phase appears with 
complex scaling properties and strong corrections to scal- 
ing, due to the presence of large internal length scales. 
This research is supported by NSF grant DMR-9700430, 
by the KOSEF through the SRC program of SNU-CTP, 
and by the Korea Research Foundation (98-015-D00090). 
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FIG. 1. Finite size scaling estimates of the critical expo- 
nents at the equilibrium point p — q for: (a) The station- 
ary state surface width W ~ L a using L — 2 L pairs for 
L — 2 4 , . . . , 2 9 . (b) The temporal surface width W ~ t 13 , 
using f/10 — t pairs at L = 2 13 . 




(a) (b) 
FIG. 2. Hill structures in the faceted growth phase, (a) 

Schematic structure with two pinning valleys. £o is the char- 
acteristic width of the hill tops, (b) Active bonds (the filled 
circles) along a local segment of a faceted ramp. 
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FIG. 3. Active bond distribution p(x) at L = 257 for 
p = 0.505, 0.51, 0.515, and 0.52 (from top to bottom). Be- 
yond foi p( x ) decays exponentially, (b) Collapse of the data 
from (a) according to Eq. (1), with /3 P = and v = 1, demon- 
strating the hills preserve their shape. 
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FIG. 4. Phase diagram of the generalized model: The 
rough equilibrium point E at r = 1 broadens into a slowly 
growing rough phase (shaded) when the digging probability r 
is reduced. The dotted line to the left of the DI critical point 
at r = represents the smooth phase. 





FIG. 5. Scaling at the p = q transition line: (a) Active 
bond distribution p(x), along p — q for system size L = 1025 
at r = 0.9,0.8,0.6, and 0.2 (ordered from the top down). 
p(x) ~ x~ x decays algebraically beyond a non-diverging hill 
top flatness length scale fo- (b) Scaling of the total active 
bond density p ~ InL/L at p = q and r = 0. 
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